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Abstract 



c^ ! We study the low-energy effective action Sefflf^] for the one-component real scalar 

field theory in three Euclidean dimensions in the symmetric phase, concentrating on 
its static part — effective potential Veff{^)- It characterizes the approach to the phase 
T^ • transition in all systems that belong to the 3d Ising universality class. We compute it 

from the probability distributions of the average magnetization in the 3d Ising model 
^ ! in a homogeneous external field, obtained by Monte Carlo. We find that the y?® term in 

H i K// is important, while the higher terms can be neglected within our statistical errors. 

Thus we obtain the approximate effective action 

^eff = d^xl -d^(pdf,(p + -m'^ip'^ + mg^Lp^ + g^ip^ \ , 

with arbitrary mass m that sets the scale, and dimensionless couplings g^ = 0.97 ±0.02 
and gQ = 2.05 ± 0.15. The value of g^ is consistent with the renormalization group 
fixed point coupling. This K//, when used instead of the traditional aip"^ + feyj^, turns 
the Ginzburg-Landau description of the long-wave properties of the 3d theory near 
criticality into quantitatively accurate. It is also relevant to the theory of cosmological 
phase transitions. 
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1 Introduction 

This work is devoted to the following problem: what is the effective potential, and the 
corresponding effective Ginzburg-Landau theory, that would provide not exact, but 
reasonably phenomenologically accurate description of the properties of the 3d Ising 
model near the phase transition (and other models that belong to the same universality 
class)? 

The model from this universality class that is particularly suitable for field-theor- 
etical treatment is the theory of one-component real scalar field in three Euclidean 
dimensions ( "3d 0^ theory" ) , defined by the (bare) action 

S = j (fx |^9^09^0 + ^mV' + A0^| . (1) 

Thus, from the field-theoretical point of view, we study the low-energy effective action 
of this theory. 

This problem, being interesting by itself, is also relevant to the theory of cosmolog- 
ical phase transitions in the early Universe. The second order high-temperature phase 
transition in the 3+1-dimensional quantum field theory is in the universality class of 
the 3d Euclidean phase transition. The weak first-order high-temperature transitions 
can be studied in the framework of effective 3d Euclidean theory as well. The effective 
potential for such problems has been a subject of recent investigations |l], §]. The use 
of the perturbation theory is hindered in three dimensions by infrared divergences and 
by the strong-coupling nature of the problem. Such issues as the existence and role of 
the |(y9p term in the effective potential remain to be settled. 

Thus the nonperturbative study of the effective action of the simplest 3d field theory 
(0), or that of the 3d Ising model, seems appropriate. 

2 The model 

We study the Ising model with the nearest-neighbour interaction on a simple cubic 
lattice. The partition function is 

Z=5]exp{/3 5]0i0, + J^0,}, 0. = ±1, (2) 

where J is the homogeneous external field ("magnetic field"). We study the symmetric 
(paramagnetic) phase, thus the coupling /3 is less than, but close to the critical value 
/3c ^ 0.22165. 

Our main subject are the long-wave (low-momentum, low-energy) properties of the 
model, when it is in the scaling region, but not exactly at the critical point. Then 
the properties are fixed, the only free parameter is the mass ( = the scale). The 
particles of the corresponding 2+1-dimensional field theory are massive (and thus can 
be nonrelativistic) and have well-defined low-energy properties, such as nonrelativistic 
scattering amplitudes. The effective action we are looking for is a convenient formalism 
to describe these properties. 
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3 The effective action 

The low-energy Ginzburg-Landau- Wilson effective action can be written as 

5e// = Jd'x i^'^Z-%cpd,ip + K//(^) - JixMx)] , (3) 

where <f{x) is the (slowly varying) average magnetization, and we keep only the lowest- 
order derivative term. To compute Seff one needs to know the effective potential 
Veff{ip) and the field renormalization factor Z^. To compute the former, it is sufficient 
to consider only the homogeneous external field J{x) = J; the latter can be derived 
from the two-point correlation function of (f. Thus we never have to work with the 
explicitly x-dependent J(x). 

3.1 The effective potential 

The effective potential (free energy) Veff{(p) [§1-|§[ is defined by 

dVeff{<^) 



dip 



J, ^= (0)j. (4) 



Thus (P), considered at the tree level, reproduces correctly the average magnetization 
(0)j as a function of external field J. We recall also that Veff{<f) is a generating 
function for one-particle irreducible (IPI) n-point Green functions r„ at momentum 
zero, 

V,ff{^)=J2-^Tr.ip = 0)cp\ (5) 

3.2 How to compute Veff on the lattice - existing methods 

Our goal is to compute Veff by Monte Carlo with the best possible precision. There 
are at least three approaches. 

1) Direct computation of the correlation functions F^, such as 

T,{p = 0) = (0^) - 3(02)^ (6) 

and so on 0-[|l. This approach works reasonably well for r4, but is hardly feasible for 
Fg and higher F„: after subtraction of disconnected and one-particle reducible parts 
the signal to noise ratio turns out to be very bad, and statistical errors are prohibitively 
high i. 

2) Compute the magnetization per spin (f = {(f))j as a function of the external field 
J by Monte Carlo, invert this function to obtain J = J{f) and integrate it numerically 



to obtain Veff{ip) according to @ 0, |T0[- The drawback of this approach is that it 
requires many measurements of (0)j at different values of J. 

3) Study the probability distribution of the order parameter -^ J2i 4>i IHl- O^is uses 



the Monte Carlo algorithm to generate a Boltzmann ensemble of configurations. For 



every configuration in the ensemble one computes the order parameter (magnetization 
per site) '^ = jfJ2i 0j; where N is the total number of sites on the lattice. Thus 
the probability distribution of if is obtained, which is characterized by the probability 
density 

To prove that this is indeed the probability density, it is sufficient to observe that for 
any function /(</?) 

Jp{^)f{v)d^={f{l^T.<t>.))- (8) 

This probability is considered usually in connection with the "constrained effective 
potential" Vconiv) M,W^-^, which is defined by 



P{^)cxexp{-QVconiv)}, (9) 

for a system in a box of volume ^2 (we set the lattice spacing equal to unity, so 57 = A^). 
It can be shown that in the infinite volume limit the constrained effective potential 
coincides with the standard Vg// |T^, [TBI . 



One can gain additional insight into these different ways to extract Veff from the 
Monte Carlo data by considering the actual distributions of the order parameter in a 
typical case (Fig. |I|). 

The first method studies, essentially, the deviation of the J = curve from the 
Gaussian, which is quite small and difficult to measure. The distribution is located at 
small values of ip, where the dominating term in Veff is (f"^, while higher terms provide 
only small corrections. 

The second method means that we use each curve to compute the expectation 
value of ip, i.e. {4>)j, and discard all remaining information contained in the data. 
(It seems that one can obtain a good method by combining this with the reweighting 
technique [|l^, |T^. Then ((/)) j could be computed accurately for any J from to ~ 0.3 



just from the data at Fig. |1| plus the data on the energy. And all the information from 
probability distributions is utilized. However, we do not pursue this line further here). 
The third method also looks at the form of the J = probability distribution and 
thus has the similar drawbacks as the first one. 

3.3 Our method 

For the Monte Carlo computation of the effective potential we have developed a method 
that is close in spirit to the constrained effective potential approach, but contains two 
significant improvements. 

First, it is obvious from the preceding discussion, that to obtain the information 
on the higher powers of ip in K//, one should study the system at ip far from zero. 
This can be achieved with the help of the external field. The definition (^ allows a 
straightforward generalization for nonzero external field: 

P{p) oc exp {-nVeoni^) + ^M ■ (10) 



Thus one can check whether some ansatz, such as Vconi'^) = fip^ + wy?^, gives a good 
approximation for the Kon, and find the values of parameters (such as r and m), by 
performing a simultaneous fit of several histograms, corresponding to different values 
of J, by (PDD, using the parameters of Vcon as fit parameters. 

The second point concerns the correspondence between the probability distribution 
and the effective potential. When we use the relation ([T^), we get the finite- volume 
constrained effective potential. However, the quantity that enters the low-energy effec- 
tive action is the infinite-volume Veff, i.e., the free energy per unit volume in the usual 
thermodynamical sense. As in the infinite-volume limit {Q -^ oo) both potentials are 
known to coincide, one is tempted to use the formula [|lj 



P{^) oc exp {-nVeffi^^) + ^Jv} ■ (11) 

and to hope that for the reasonably large volume Q the finite volume effects are reason- 
ably small. It turns out that this is actually not the case. The situation is illustrated by 
the lower graph at the Fig. ^. The finite volume corrections (which manifest themselves 
by the volume dependence of parameters of K//) are unacceptably large, in spite of 
going to zero in the infinite volume limit. This behaviour shows that (O) is not a very 



good approximation. It turns out that while (|TI]) reproduces correctly the exponential 
dependence on ^2 at fi — > oo, it misses the important preexponential factor. 

The improved formula for the probability distribution of the order parameter of the 
system in a finite box of volume Q with the periodic boundary conditions is 



Pi^) « \j^-^^^ exp {-nv^ffi^) + nM , (12) 

where Veff {(f) is the infinite- volume effective potential. This relation can be found 
(in somewhat implicit form or for special cases) in the literature (see \W\; |T^,(5.54); 



|T^,(4.27)). The upper graph at Fig. ^ confirms that this is indeed a good approxima- 
tion: no finite volume corrections are visible for L/^ > 4. 

We outline here briefiy the derivation of the preexponential factor in (|1^). On the 
one-loop level one can write the effective potential as a bare potential plus one-loop 
correction. For a theory with a bare Lagrangian 
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Co = 7;d,(Pd,^ + Vo{(P) (13) 



we have 

K„(.)-f»(.)4/(gj.n(p= + 0). (14) 

At the same time for the constrained effective potential in the box of volume fl the 
integral is substituted by a sum over momenta allowed by the boundary conditions, 
with the p = mode excluded: 



Thus the difference between Vcon and Vg// reduces to the difference between the sum 



and the integral. This difference has been computed in []T8[, sect. 3.3. The leading 
large volume term takes a simple form: 



— ln^° 
2fi "^ difi'^ 



Vcon{^)\n-Veff = -7^iri--^. (16) 



This provides the preexponential factor in (|T2|). The last step is to substitute Vq by 
Vgff, which means essentially the use of selfconsistent values for the masses. 

3.4 The preexponent in the probabiUty distribution 
and the equal weight versus equal height problem 
for asymmetric first order phase transitions 

We deviate a little from our problem to make an additional comment on the for- 
mula (0). 

This formula seems quite universal and useful for the various problems connected 
with the order parameter probability distribution. In our problem we use it in the 
symmetric phase for a range of J. A very different and in some sense complementary 
application is the equal weight versus equal height problem for asymmetric first order 
transitions Il3)[01"lHI- Consider a model with a first order phase transition between 



two phases characterized by the order parameters (pi and (p2, exactly at the transition 
point. This means that the free energies for both phases are the same, i.e., Veff{fi) = 
Veff{^2)- Put the system in a finite box much larger than the correlation lengths in 
both phases. Consider the probability distribution P{<f) of the order parameter. Then 
this distribution can be approximated by the two narrow Gaussian peaks around ipi 
and {p2 pO|, [l^. The question is, do they have equal height or equal weight? The 
formula (pTf ) predicts equal height, while it can be shown that equal weight is the 
correct answer pl|-||23[,[0. But this it just what follows from (p!2D. 



4 Monte Carlo computation of the effective action 

Now we turn to our computation of the effective action. We study the 3d Ising model (H) 
on a simple cubic lattice with periodic boundary conditions, on lattices from 14^ to 
58^. The Swendsen-Wang cluster Monte Carlo algorithm in the external magnetic 
field [Q-[^ is used to generate the Boltzmann ensemble of configurations. (We use the 
version of this algorithm without the ghost spin). For every configuration we measure 
magnetization per site ip = j^J2i4>i and compute the histograms for the probability 
density P{<f), for several values of J. Then we do the simultaneous fit of all the 
histograms with the formula (|12D- (We minimize the sum of x^ from the individual 
histograms). The first ansatz to try is 

Veffi^) = V + W\ (17) 



inspired by the standard Ginzburg-Landau theory (or the tree-level 0^ theory). The 
result is shown at Fig. |^. One can see the discrepancy between the data and the fit 
— not very large, but statistically significant at our level of precision. This means 
that (|r^) does not provide a good quantitative description of the effective potential. 
So we consider a three-parameter expression 

Ve/zM = V + V + W;</'^- (18) 

This form is motivated by several intuitive considerations that are discussed below. We 
have found that it provides the ideal fit (Fig. |I|): there is no systematic discrepancy 
between the data and the fit, just statistical noise. We have found no other reasonable 
ansatz that works so well (we have also tried three-parameter expressions V"e//(v?) = 
np^ + mIv^I"' and Veff^f) = np"^ + w\ipf + mp^ — they work poorly). 

Thus for every value of the bare coupling (3 we obtain the low-energy effective 

Lagrangian 

1 

2' 

The three parameters r, m, w are determined by the fitting procedure described above. 
The field renormalization factor Z^p is obtained from the propagator 



£e// = -Z-'d^^d^^ + r^^ + u^^ + w^\ (19) 
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We use the lattice 


version of this. 










G2(P)- 


-' = Z-. 


^>;2(i- 
11=1 


COSPfj 


) + 2r 



(20) 
(21) 

(22) 
(23) 



and use typically three values of momentuin, 



'2,7m 
p^=(p,0,0), P=-^, n = 0,l,2, (24) 

to find Z~^ as the slope of G2(p)~^ as a function of J2fi 2(1 — cosp^). 
After the renormalization of ip, 

V = \JZpVRj (25) 

we obtain the effective Lagrangian in the form 

^eff = i^df.^ndf.^R + -mV?j + mg^^^ + gQ^% (26) 



where 

m = j2Z~r, ^4 = ^ — , Qq = Z^w. (27) 

^ y2V- 

In the continuum hmit (m ^ 0) this effective Lagrangian should be universal. Thus 
the only free parameter is m, that determines the scale, while the dimensionless four- 
and six-point couplings g^ and qq take definite values that are the same for the whole 3d 
Ising universality class. Our numerical results are collected in Table |I| and represented 
at Figs. |-0. 

5 Data analysis and extrapolation 
to the continuum limit 

Apart from statistical errors, there are two sources of systematic errors: finite volume 
and finite UV cutoff. To check for the finite volume effects, we increase the lattice size 
L, keeping m fixed. One can see from Fig. |^ that for L/^ > 4 finite volume effects are 
negligible. To check for the effect of the finite UV cutoff, we keep L/C, fixed at ^4.1, 
increase ^ and scale J according to 

J oc C^^^'^ iPS ^ 1.57, z/ ^ 0.63) (28) 

The scaling limit is characterized by stabilization of g^ and qq as the functions of ^, 
for large enough ^. We observe a smooth approach to scaling (Figs. § and ^, but note 
that the shift of Qq caused by the finiteness of the cutoff is still visible on our largest 
lattices. The correlation length in zero field is ^ = m~^ ^ 14 for the 58^ lattice, but 
one should keep in mind that it is smaller for nonzero J, and for the largest J used in 
the fit it is approximately twice as small. The presence of finite cutoff effects in our 
data makes it necessary to extrapolate the results to continuum limit ^ — i> oo. As we 
keep L/C, fixed, we work in terms of L. The reasonable extrapolation is 

geiL) = gei^x^) + aL-\ (29) 

The statistical errors in the data for gg are still too high to allow the determination 
of the exponent k with reasonable accuracy, while different values of k (such as k = 1 
or K = 2) lead to considerably different extrapolated values gQ{oo). The situation can 
be alleviated if we take into account that the statistical errors of g^ and g^ are not 
independent (Fig. ^, and so there is a linear combination of g^ and g^ that has much 
smaller error than (74 and gQ separately (Fig. [Zp. For this combination one can find 
the corresponding exponent k = 1.5 ± 0.2. The plausible assumption is that g^ should 
be extrapolated with the same exponent {g^ shows little dependence on L for L > 17, 
so one does not really have to extrapolate it). That is why we plot the couplings at 
Figs. ^ and ^ as functions of L~^'^. 

Our result for the continuum limit is 

^4 = 0.97 ±0.02, 

c/6 = 2.05 ±0.15, ^^ 

where the errors are the standard deviations. 



6 Discussion 

Here we compare our results with the data available in the literature, and give some 
semiintuitive arguments in favour of the effective Lagrangian (^6]) as a good approxi- 
mate theory. 

6.1 Four-point coupling 

The four-point coupling of the 3d scalar field theory has been a subject of many studies. 
It is not the main subject of our investigation (our main point is the role of Qq), but 
the comparison of our (74 with the available data provides a useful consistency check of 
our computation. Here we list some data for (74 obtained by different methods. 

From the 3d renormalization group fixed point: 
0.989 ±0.004 [g3]-[§U| 



0.981 [plU, revised estimate 



From the high-temperature series ||32||-pl^ 

0.99 ± 0.03 (simple cubic lattice) 
0.991 ± 0.003 (BCC lattice) 
0.99 ± 0.01 (FCC lattice) 

From the Monte Carlo studies of the (0^) — 3(0^)^: 
1.00 ±0.04 11,^ = 3.3 
1.00 ±0.03 0,^ = 6.6 
0.80 ±0.02 0,^ = 16 
0.95 ±0.08 §,^ = 3.2 
0.90 ± 0.04 p5|, extrapolation to the continuum limit 



1.00 ±0.08 111,^ = 14.5 



We observe that our result g^ = 0.97 ± 0.02 fits well into this picture. 

6.2 Six-point coupling 

The information about qq is more scarce. The only Monte Carlo study we are aware of 
is |§]. However, large statistical errors made it impossible to reach a definite conclusion 
about the value of Qq in the continuum limit and whether it is different from zero. 

Another source of information is provided by the study of the "Ising equation of 
state" in the framework of the e-expansion [|5^-|^. The equation of state describes 



magnetization as a function of the homogeneous external field, thus providing infor- 
mation on the effective potential. As the gradient term in Ceff, that determines the 
normalization of the renormalized field, is not considered, one can extract parameters 
that do not depend on this normalization. The ratio g^/gl is such a parameter invari- 
ant under the change of the scale of the field. The convenient representation of the 



equation of state is given by Avdeeva and Migdal [37|. In original notation 



^ = ^-(/34-7)/70(Mx^/7), (31) 



where M is magnetization, H is the external field, x = dM/dH is the susceptibility in 
the finite field, and the function 6 is found to be 



(f){m) = m — m^ H m^ + 0{e^). 



(32) 



for the dimension of space d = 4 — e. From this equation it is straightforward to 
compute 

96 . 20 



W 2.-^.^ + 0(.3). 



(33) 



The next order in e can be derived from the parametric representation of the equation 



of state described in [^], sect. 25.1.3. One finds 

20 



96 



(9.) 






(34) 



The series is obviously divergent, and the only straightforward conclusion for e = 1 
that one can make without some resummation is 



96 

i9.y 



2.0 ±0.7. 



(35) 



This agrees well with our result (|50|) . 

The effective potential of the 0'^ model can be also computed in the framework of 
perturbation theory directly in d = 3. The result including up to five loops can be 
found in ^^, eq.(3.3). It turns out, however, that the two-loop contribution to Qq 
is four times as large as the one-loop, the three-loop even larger, etc. This makes it 
impossible to derive any number for qq without some resummation of the series. 

Yet another approach to the computation of the effective Lagrangian is provided 
by the Wegner-Houghton equation in the local-potential approximation [^, |^. The 
fixed point for the potential corresponds to qq = 2.40 p4 |. 

One more approach is the strong-coupling expansion [0-[0 and the dimensional 
expansion |^ for the field theory in the Ising limit, i.e., in the limit of infinitely strong 
bare coupling. The results obtained ioi d = 3 disagree with ours, as the strong-coupling 
expansion favours 9q = P7|, while the dimensional expansion leads to conjecture that 



oo 



Finally, we compare our qq with results of Tetradis and Wetterich, obtained within 
the "effective average action" approach [^-|51|. They use the symbol u^ for the six- 
point coupling, the correspondence is g^ = u^/iS. The fixed-point value is M3* = 
87.4 (|§0[' Table 2), while the asymptotic value of the low-energy coupling, as the 
phase transition is approached from the symmetric phase (which is more appropriate 
to compare with our g^), is Uss = 107 |S^. This corresponds, respectively, to g^ = 1.82 
and ge = 2.23. 



6.3 Philosophy 

A widespread point of view on the effective potential in 3d is as follows. The problem 
should be considered in the framework of the 0^ theory. Then either you work on the 
tree level, and have a standard Landau theory with 

V,ff{y,) = ay,' + by,\ (36) 

or you include loop corrections, and then you must retain all powers of ip in Ve//. The 
ip^ term should be considered on equal footing with other higher terms. 

Our study corroborates an alternative point of view advocated by Tetradis and 
Wetterich |^0|, ^: that while (^) is a rather rough approximation, the ansatz 

"^effiv) = 2"^V^ + mg^Lp'^ + geip^ (37) 

gives a very good approximation, and the higher powers of ip can be considered as 
small corrections. Additional arguments in favour of this point of view can be drawn 
from the smallness of the critical index rj in the 3d theory. 

Consider the asymptotic behaviour of magnetization (p = {(f))j in the external field 
as a function of J at large (p ( = large J). This is the same limit as if we keep J fixed 
and let m — »• 0. The equation dVeff/dip = J gives 

Gge^" = J, {<l))j oc J'/^ (38) 

The Landau theory ( p6l) at the critical point (a = 0) would give 

ip^ oc J, {(j))j oc J^/l (39) 

However, this is exactly the definition of the critical index 5: 

{<P)j ex J'l'. (40) 

Thus the Landau theory ( p6D corresponds to 5 = 3, and (^) corresponds to 5 = 5, 
while the correct value is 

5 =^±^^-4.8. (41) 

We see that while 5 7^ 5 means that ( p7| ) is, strictly speaking, inapplicable for very 
large y? (such that V'J^^^{}p)) ^ K//(0) — ^), 5 is still a much better approximation to 
S than 3, due to the smallness of rj. In the large N limit of the 0{N) Heisenberg model 
in 3d ?7 ^ 0, and the A^-component analog of ( ^7|) becomes exact [^ ^ . 



The form {p7\ ) is inapplicable for m = 0, i.e. exactly at the transition point. In this 
case every nonzero ip is 'very large', and the correct form is 

Veff(Xip'+' (m = 0). (42) 

However, the whole idea of the low-energy effective action is no more applicable in this 
case, because 'low energy' means p -C m. So we always work at nonzero m. 

Notice that cp^ is renormalizable in 3d, while higher terms are nonrenormalizable. 
Thus ( ^71) amounts to the idea to keep just all renormalizable terms in the effective 
action. The same phenomenon — that this is a good approximation — is observed for 



the 3d 3-state Potts model [^, and seems to be rather general. 
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7 Conclusions 

We have studied the low-energy effective action for the scalar field theory in three 
dimensions. We considered the 3d Ising model in the symmetric phase in the scaling 
region. We have found that a very good approximation is provided by 

Veffi^) = -m^if^ + mg^ip^ + geif'^, (44) 

and estimated the values of dimensionless couplings: 

^4 = 0.97 ±0.02, ^6 = 2.05 ±0.15. (45) 

(The errors are standard deviations). This effective action is universal for the whole 
3d Ising universality class, in the sense that the only free parameter is the mass m. Its 
applicability region is bounded by the requirement that cp is not too large. It works 
well at least for cp such that meff{p) = [Veffif)]^^"^ — 3"^, and maybe further. 

We observe a smooth approach to scaling in the 3d Ising model. Our results for 
g^ and g^ are in good agreement with the values available in the literature, but dis- 



agree with [^-[^. While, to our knowledge, our result for g^ is the most precise 
among the available Monte Carlo estimates, and this is the first Monte Carlo study to 
provide the definite result on gg, we consider the direct Monte Carlo check that the 



approximation (|^) works perfectly as our main result. 

This is in agreement with the results of Tetradis and Wetterich, obtained within 
the effective average action approach fSOl, pll], that (p^ is important, while the higher 



terms in Vejf provide only small corrections. 

Our main tool was the study of the probability distribution of the order parameter 
in the finite system with the periodic boundary conditions, in the external field. We 
find that the simple asymptotic formula ([T^) that connects the probability distribution 
in a large, but finite system with the infinite-volume effective potential works perfectly 
in our case, when L/^ > 4. The formula seems to be rather general and useful also for 
other problems in the theory of the order parameter distribution. 

The Ising model in the symmetric phase is just the simplest three-dimensional 
model. It would be interesting to extend the approach developed here to its broken 
phase, to weak first order transitions and to models that serve as effective 3d theories 
for high-temperature phase transitions in QCD and in the Higgs sector of the standard 
model. 
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L^ 


14a 


17^ 


22^ 


30^ 


38^ 


58^ 


P 


0.2115 


0.2142 


0.2167 


0.2186 


0.2195 


0.22055 


J 


0.0 


0.0 


0.0 


0.0 


0.0 


0.0 




0.0043 


0.0027 


0.0014 


0.00066 


0.00036 


0.00013 




0.013 


0.0081 


0.0043 


0.002 


0.0011 


0.00038 


Nconi 


3 X 700000 


3 X 720000 


3 X 700000 


3 X 400000 


3 X 500000 


3 X 100000 


y-l 

^(J=0) 


0.2294(2) 


0.2351(2) 


0.2407(2) 


0.2458(3) 


0.2493(3) 


0.2550(6) 


r 


0.01014(2) 


0.00690(3) 


0.00414(1) 


0.002263(6) 


0.001468(5) 


0.000639(7) 


u 


0.0163(2) 


0.01337(20) 


0.01062(11) 


0.00814(14) 


0.00657(12) 


0.00451(17) 


w 


0.0394(7) 


0.0388(8) 


0.0381(5) 


0.0360(10) 


0.0361(12) 


0.0343(20) 


Ex' 


188 


208 


193 


215 


200 


140 


-^^bins 


165 


161 


161 


229 


160 


146 


m 


0.2973(3) 


0.2423(5) 


0.1855(3) 


0.1357(2) 


0.1085(2) 


0.0707(4) 


94 


1.041(12) 


0.998(15) 


0.988(10) 


0.993(17) 


0.974(18) 


0.98(4) 


56 


3.26(6) 


2.99(6) 


2.73(4) 


2.42(7) 


2.33(8) 


2.07(12) 



Table 1: The numerical results obtained by Monte Carlo. A'^conf is the number of 
configurations used. J^x' is the sum of x^ for all histograms (minimized by the fit); 
A^bins is the total number of bins in the histograms. The numbers in the parantheses 
are standard deviations of the last decimal digits. 
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Figure 1: A set of probability densities P{f) for the magnetization per lattice site cp, 
for the Ising model (0). The solid line is the fit with (|T2p, Veff{(p) 
Three histograms are fitted simultaneously with the same Vgff. 
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Figure 2: Finite volume dependence of the parameter r, obtained by fitting the 
probabihty distributions for f3 = 0.2186, J = 0, 0.00066 and 0.002 — the same values 
as at Fig. |T] — with (|TTp and ([T^), Veff{(p) = np^ + uip^ + wip^, for different lattice 
sizes: 20^, 24^, 30^ and 38^ (the correlation length ^ 
less volume-sensitive. 
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Figure 3: The data is the same as at Fig. ^ The solid line is the best simultaneous fit 
with (p!2D, Veff{<~p) = TLf^ + u(f'^. Statistics: 400000 configurations for every histogram. 
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Figure 4: The dimensionless four-point coupling g^ as a function of the lattice size L. 
The ratio L/C, is kept fixed at around 4.1. The errors shown are standard deviations. 
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Figure 5: The six-point coupling qq as a function of lattice size L. {L/^ ^4.1^ 
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Figure 6: The simultaneous plot of (74 and qq. The ellipses show the statistical 
errors (semiaxis = standard deviation). The lattice size L is shown near the ellipses, 
L/^ Ri 4.1. The broken line is for the same parameters {13 and three values of J) as 
for L = 30, but measured on a 38^ lattice. This provides the additional check for the 
smallness of the finite volume corrections. 
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Figure 7: The linear combination of dimensionless couplings that is determined with 
the smallest statistical error. 



